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Abstract. Inflationary homogeneous isotropic cosmological models filled by scalar 
fields and ultrarelativistic matter are examined in the framework of gauge theories of 
gravitation. By using quadratic scalar field potential numerical analysis of flat, open 
and closed models is curried out. Properties of cosmological models are investigated in 
dependence on indefinite parameter of cosmological equations and initial conditions at 
a bounce. Fulfilled analysis demonstrates regular character of all cosmological models. 
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1. Introduction 

As it was shown in a number of papers (see m and Refs, given therein), general feature 
of homogeneous isotropic cosmological models built in the framework of gauge theories of 
gravitation (GTG) [Hill] is their regular bouncing character. All cosmological solutions 
of generalized cosmological Friedmann equations (GCFE) deduced for homogeneous 
isotropic models in GTG are regular in metrics, Hubble parameter and its time derivative 
because of gravitational repulsion effect at extreme conditions (extremely high energy 
densities and pressures), which takes place in the case of usual gravitating systems 
satisfying energy dominance condition j[H|. In the case of inflationary cosmological 
models filled by interacting scalar field and usual gravitating matter, indicated above 
property of cosmological solutions was obtained without concretization explicit form of 
scalar held potential and equation of state of gravitating matter Tj. However, we have 
to use their explicit form to obtain and analyze cosmological solutions of GCFE. 

The present paper is devoted to analysis of inhationary cosmological models filled 
by noninteracting scalar helds and ultrarelativistic matter. In Sec. 2 the GCFE for such 
inhationary cosmological models are introduced. In Sec. 3 some general properties of 
its solutions are examined. In Sec. 4 and 5 by using quadratic scalar held potential 
numerical analysis of cosmological solutions for hat, open and closed models is carried 
out. 
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2. Generalized cosmological Friedmann equations in GTG 
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The GCFE deduced in the frame of Poincare GTG and metric-affine GTG have the 
following form [6-8] 
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where R(t) is the scale factor of Robertson-Walker metrics, p is energy density, p 
is pressure, k = +1,0,—1 for closed, flat, open models respectively, G is Newton’s 
gravitational constant (the light velocity c = 1), a is indefinite parameter with inverse 
dimension of energy density, a dot denotes differentiation with respect to time t. 
Together with Newton’s gravitational constant G the parameter a is fundamental 
physical constant of relativistic isotropic cosmology built in the frame of GTG. The 
value of C+ 1 determines the scale of extremely high energy densities [SJ, which will be 
supposed to be less than the Planckian one. From equations (l)-(2) follows the energy 
conservation law in usual form 


p + 3H (p + p) = 0 



(3) 


Obviously, the solutions of GCFE m-m depend on values of G and a. However, we 
can exclude them and transform equations 0-0 to dimensionless form by means of 
the following transformation of the time t and functions R(t), pit), p(t) 
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where a prime denotes the differentiation with respect to t. The form of equation m 
does not change by transformations 0 


p + 3H (p + p) = 0 



(7) 


In the case of models filled by scalar field 0 with potential V 
matter with energy density p m and pressure p m we have 

P = ^ 0 2 + V + p m , p = ^ 0 2 - V + Pm- 


V (0) and gravitating 
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By using the transformation © we obtain from (JHJ) the following dimensionless 
expressions: 

P = 2 + V + /5 m , p = - V + p m , (9) 

where <p = (pVG, V = aV, p m = a p rn , p m = ap m . By neglecting the interaction 
between scalar held and gravitating matter, from conservation law 0 follows the 
equation for scalar held 


0 " + 3 = 

o<p 

and conservation law for gravitating matter 


( 10 ) 


Pm + 3 H(p m + Pm) ~ 0. 


( 11 ) 


To integrate equations (5)-(6) with p and p defined by (9) we have to know the 
equation of state of gravitating matter p m = p m (p m ). To build realistic cosmological 
models we have to take into account the change of equation of state under their evolution 
in accordance with properties of gravitating matter, but it is generally a complicated 
physical problem. To study cosmological models near a bounce and simplify further 
analysis we will use below the equation of state for ultrarelativistic matter p m = \p m . 

In considered case from equation (11) follows the integral for gravitating matter 

pR 4 = const. (12) 

Then the GCFE (5)-(6) can be transformed in accordance with Ref. [5] to the following 
form: 
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where Z = 1 + 4V — f' 2 . 

The greatest part of solutions of GCFE (13)-(14) are inhationary cosmological 
solutions having direct physical interest. Properties of these solutions are studied below 
in sections 3-5 in dependence on initial conditions at a bounce, restrictions on parameter 
a and parameters of scalar held potential. Note that, if the value of a -1 is much less than 
the Planckian energy density, equations (13)-(14) with scalar held potentials applying 
in inflationary cosmology lead also to oscillating solutions studied in Ref. HI • 
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3. Some properties of inflationary cosmological solutions 
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Let us introduce 4-dimensional phase space P with the axes (0, cfl , tilH , p m ). Then some 
trajectory in space P corresponds to any solution of GCFE. Trajectories of cosmological 
solutions are situated in area of P, where Z > 0 and, hence, are limited by the bounds 
L defined as 

Z = 0 or 0' = ±Vl + 4H. (15) 


According to equation (13) values of the Hubble parameter are equal to 
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Two values of the Hubble parameter H + and 77_ correspond to two different solutions 
named as 77±-solutions. In points of bounds L we have D = 0, 77 + = IP and the 
Hubble parameter is determined by 

_1_ 9V_ 

0' (90 

Now let us consider extremal T7o-surfaces in the space P, in points of which 77+ = 0 
or 77_ = 0. Denoting the values of quantities on 77 0 -surfaces by means of index “0”, 
we obtain from (13) the following equation for T7o-surfaces, which are in fact in 3- 
dimensional hyperspace 77 = 0 of the phase space P: 
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where Z 0 — 1 + 4Vj — 0 q 2 . For given scalar field potential V (0) this equation determines 
77 0 -surfaces for flat models (k — 0), and in the case of closed and open models equation 
(19) gives parametric families of 77 0 -surfaees with parameter R 0 . According to equation 
(14) the time derivative of the Hubble parameter 77g on 77 0 -surfaces is 
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If 77' > 0, we have in corresponding points of 77o-surfaces the transition from 
compression to expansion (a bounce), and in the case 77 q < 0 the transition from 
expansion to compression takes place. Equations (19)-(20) permit to End specific points 
of solutions, where 77 = 0 and 77' = 0 simultaneously. 
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To analyze solutions of GCFE we have to study H 0 -sur faces defined by (19) for 
given scalar field potential. We will use below the simplest quadratic potential in the 
form V = \7n 2 qi 2 or V = where fh = In this case H (r surfaces depend on 

parameter fh and the value of scale factor R 0 . In the case of p m = 0 equation (19) for 
i/o-surfaces includes only variables </>o and ft 0l and the analysis of //o-surfaces as well as 
bounds L is reduced to analysis of corresponding curves on the plane (</>, Numerical 
analysis will be given below in the Planckian system of units with h — c — G — 1. Note, 
that admissible values of dimensionless potential V = aV can be essentially greater 
than 1, if the scale of extremely high energy densities defined by or 1 is much less then 
the Planckian one. Hence, values of parameter fh can be large. Below the analysis is 
given for various types of cosmological inflationary models. 

4. Numerical analysis of flat inflationary models 

In the case of flat models (k = 0) // 0 -curves B i2 are presented in Fig. 1 together with 
bounds L\ 2 having the form of two hyperbolas ft = ±\J 1 + 2 fh 2 f 2 . Each of two curves 
B\ 2 contains two parts corresponding to vanishing of H + or and denoting by (Bi + , 
B 2+ ) and (Bi_, B 2 -) respectively [I], The parts (_B 1+ , Bi_) and (B 2+J B 2 _) have the 
common points (0 = 0 ,ft — ±1) with the bounds L\ and L 2 on the axis ft. From 
equations (19) - (20) it is possible to show, that if 

fh < 2^-3 - 47 (l009 + 84\/l59) + (l009 + 84\/l59) 1/3 « 1.346, 

than the value of H' Q on i/o-curves is positive (this means that H' +0 > 0 in points of B\ + 
and B 2+ , and H'_ 0 > 0 in points of Bi_ and B 2 _). In this case ^-curves are “bounce 
curves” [I], in points of which the transition from compression to expansion takes place. 
If m > 1.346 (that corresponds to large values of a), the regions with Hq < 0 appear on 
// 0 -curves near the axis ft for sufficiently small values of \<j>\ (see Fig. CD. In this case, 
besides regular cosmological solutions, the GCFE have scalar field oscillating solutions 

El- 

Now let us analyze properties of cosmological solutions for flat models. Each such 
solution, generally speaking, includes the quasi-de-Sitter compression and inflationary 
stages, stage after inflation (Fig. [2D and the stage of transition from compression to 
expansion (Fig. [2D- Properties of cosmological solutions — duration of different stages 
of evolution, characteristics of scalar field oscillations and behaviour of the Hubble 
parameter after inflation and before the compression stage are studied in dependence 
on initial conditions at a bounce and the value of fh. Like to GR na, the initial value 
of cf )o must be greater than 1 (0 O > 1) to appear inflationary stage. 

Because GCFE lead to slow-rolling regime which coincide practically with that of 
GR, properties of inflationary stage of our models are the same as in GR. Namely, the 
duration of inflationary stage t in fl increases linearly by increasing of \<f 0 \ and decreases 
inversely proportional by increasing of fh (Fig. 0. 
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Figure 1. Bound Ti, 2 -curves and bounce curves B 12 for flat models in the case of 
potential V = ^ifi 2 f 2 . Dotted curves correspond to negative values of H' 0 . 
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Figure 2. Compression and inflationary stages; postinflationary stage. 
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Figure 3. Transition stage 
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Figure 4. Duration of inflationary stage in dependence on values 0 O and to. 


Because of 0(, 0, cosmological solutions are asymmetric in the time with respect 

to the bounce point. In particular, the duration of quasi-de-Sitter stage of compression 
is not equal to duration of inflationary stage. In the case of bounce in points of curves 
B i + and B 2+ the duration of inflationary stage is greater than duration of quasi-de- 
Sitter compression stage. If the bounce takes place in points of bounce curves B\_ and 
B 2 - the situation is opposite. 

For given value of fh inflationary models without quasi-de-Sitter compression stage 
{p m = 0) can be realized for certain interval of 0o- For example, at fh = 1.2 
corresponding interval of 0 O is (28, 100), then the growth of the scale factor R{t ) is 
in interval (e 130 , e 1450 ). By increasing of fh the upper and lower bounds of interval for 
0o °f discussed solutions increase linearly (Fig. Or)- hi the case of presence of radiation 
{pm ~f~ 0) the lower bound of interval of 0 O for solutions without quasi-de-Sitter stage of 
compression decreases essentially and upper bound decreases slowly (Fig. EJo) . 

The most principal feature of our cosmological models is the presence of regular 
transition stage from compression to expansion. The behaviour of scalar held 0 and the 
Hubble parameter H at transition stage is essentially nonlinear and complicated. The 
duration of transition Af tr stage decreases by increasing of each of parameters fh, 0o 
and p m0 (Fig. El - Fig.d). 

After inflation stage the scalar held oscillates near its minimum. These oscillations 
don’t depend on initial values of 0 O and p m0 , but essentially depend on parameter 
fh. Like to GR frequency and amplitude of scalar held oscillations decreases with the 
time, but quantitatively such dependence is different. Note, that energy density p m is 
negligibly small at the end of inhationary stage, even if p m0 0. Application of Krylov- 
Bogoliubov method m of solutions of the system of differential equations (10) and (14) 
with p m = 0, after some tedious calculations, leads to the following asymptotics for 
scalar held 0(f) and the scale factor R{t) at late time 



( 21 ) 
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Figure 5. The upper and lower bounds of interval for admissible values of <j>o for 
solutions without quasi-de-Sitter compression stage. 


Attr 




Figure 6. Duration of transition stage in dependence on values of rh and <p q. 
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Equation m shows, that the amplitude of scalar held oscillations A decreases with time 
as follows A = a{rh)t~ x , where a = 4?7i _1 /a/ 3. Numerical analysis conhrms relations 
(ED-© at sufficiently small values of m {rh < 30), in particular, the amplitude 
of oscillations decreases inverse proportional to time (at rh = const) and inverse 
proportional to rh (Fig. EJ)4 But at large values of rh (rh ~ 1000) conclusions obtained 
by means of Krylov-Bogoliubov method differ from numerical results. In particular, the 
inverse proportional dependence between the amplitude of scalar held oscillations and 
the value of rh is not valid (Fig. EJ). This is connected with the fact, that in the case of 

f According to numerical calculation, constant of proportionality at F ig. |HJo is equal to 3.070, which is 
close to y/Srr « 3.069 as Krylov-Bogoliubov method requires. 
















Analysis of inflationary cosmological models in gauge theories of gravitation 


9 


At tr 



10 xp m o 


Figure 7. Duration of transition stage in dependence on p m o- 


large values of rh (i.e. large values of a) the Hubble parameter, unlike to GR, oscillates 
near the value H = 0 fFig. HOI) . The frequency of //-oscillations is two times greater 
than the frequency of 0-oscillations. The value of H becomes positive with the time 
and tends asymptotically to zero like to GR. Note that, there are similar 0-oscillations 
and //-oscillations before the compression stage also. 

~-i _-i 
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Figure 8. Dependence of amplitude of 0-oscillations on the time t and on rh at 
m < 30. 


5. Numerical analysis of open and closed models 

In the case of open models //o-curves are similar to that of flat models and the family 
of Z/ 0 -curves is situated on the plan (0,0') between curves L x and Z?i, and also between 
curves L 2 and B 2 . At Ro —■> 0 //o-curves of open models tend to bounds L i >2 , and at 
Rq —> oo they approach to curves Bi 2 (see Fig. 1). In accordance with expression (21) 
for //' the minimum value of rh , by which negative values of //(, appear for open models, 
depends on R 0 and it is smaller then that for flat models. Properties of cosmological 
solutions for open models are similar to that of flat models. There are area of H' < 0 
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Figure 9. Dependence of f- 
oscillations amplitude at large 
values of th. 
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Figure 10. iF-oscillations after 
inflation at large values of m. 


on i/o-curves, if m > 2ftn/6 ~ 1.45 and the value of Ro is arbitrary or m < 2ftn/6 
and 

~ [Y I 1 + 24?n 2 
V 47T y 1 + 30m 2 + 216m 4 

Note that in the case of open models there are also solutions with Z < 0, which are 
separated from cosmological solutions by bounds Li and L 2 . These solutions are hardly 
ever have physical sense. 

The form of FFo-curves for closed models (k = +1) differs essentially from that for 
flat and open models, and it depends on the value of R 0 for given value of m. For 
sufficiently large values of Ro H 0 -curves for closed models include besides two curves 
similar to that of flat model also one closed curve with center in origin of coordinates 
(Fig. ED. If m < 1.34 and the value of R 0 is sufficiently large, the value of H' 0 is positive 
on i/o-curves (a) and on the part of Ho-curve (6). In the case m > 1.34 the region with 
K < 0 appears on H 0 - curves (a) near the axes ft . By decreasing of R 0 the // 0 -curves 
of closed models are deformed and transformed into Ho-curves presented in Fig. El If 
m < 1.4 L/g-curves (c) contain regions with H’ Q > 0 and regions with H' 0 < 0; small 
region with H' 0 < 0 can be present on FF 0 - curves (d) although on the greatest part of H 0 - 
curves ( d ) we have Hq > 0. In the case m > 1.4 and Ro < Rqi (for example, i?oi = 1-5 
for m = 2.1, Rqi = 2.5 for m = 3.5)§ we have H' 0 > 0 on the curves (/) and H' 0 < 0 
on the FFo-curves (e). Like to flat and open models oscillating solutions exist for large 
values fh (a). Unlike to flat and open models, in the case of closed models there are 
cosmological solutions symmetric in the time with respect to a bounce ( ft 0 = 0). Note, 
that cosmological solutions have properties similar to flat models, if the value of Ro at 
a bounce is sufficiently large. 

§ The value R ^ is the real root of cubic equation, which can be obtained from equation (19) 
tt(-243to 6 - 3067nh 4 - 4087r 2 ?h 2 + 32 t t 3 )R^ + 3(6561m 8 + 7776 tt m e + 1296C>7r 2 m 4 - 5376 tt 3 m 2 + 
2567r 4 )i?g 1 — 6912m 2 7r(81m 4 — 1447 tto 2 + lGn 2 )^,^ + 3981312fh 4 7r 2 = 0 
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Figure 11. H 0 -curves for closed models at large values Rq- 
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Figure 12. Ho-curves for closed models at small values of Ro- 


6. Conclusion 

Numerical analysis of inflationary cosmological models curried out in the framework of 
gauge theories of gravitation confirms obtained earlier conclusion about their regular 
character. The difference in behaviour of cosmological models at the end of inflationary 
stage in comparison with GR (in the case of large values of a) can be interesting, 
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if we take into account recent results concerning the anisotropy of cosmic microwave 
background [EJ. This means that further investigations of perturbations in discussed 
inflationary models is of direct physical interest. 
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